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Calculation of Unsteady 
Linearized Euler Flows in 
Cascades Using Harmonically 
Deforming Grids 

Kenneth C. Hall 1 2 
William S. Clarks 

ABSTRACT A method for calculating unsteady, in viscid, compressible 
flows in cascades is presented. Using the linearized Euler technique, the 
flow is decomposed into a steady or mean flow plus a harmonically vary- 
ing small disturbance flow. The equations that describe the small distur- 
bance flow are linear variable coefficient equations, and are solved using a 
pseudo-time time marching Lax-Wendroff technique. Unlike previous lin- 
earized methods, however, the solution is computed on a harmonically de- 
forming computational grid that conforms to the motion of the vibrating 
airfoils. The mean flow and perturbation flow solutions are defined in the 
deforming coordinate system rather than in a coordinate system fixed in 
space. Hence, no extrapolation terms are required to implement the upwash 
boundary conditions at the airfoil surfaces significantly improving the accu- 
racy of the method. For transonic flow calculations, unsteady shock motions 
are modelled using shock capturing. The unsteady loads due to the shock 
motion are then seen as pressure impulses. Representative computational 
results are presented for transonic channel flows and subsonic and transonic 
cascade flows. 


1.1 Introduction 

To understand and predict flutter and forced response in turbomachinery 
blade rows, aeroelasticians require accurate and efficient models of the un- 
steady flow fields that arise from blade motion and incident gusts. Over 
the past two decades, many computational methods have been developed 
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to analyze unsteady flows in turbomachinery. For example, the linearized 
potential methods have been developed extensively [WG81, VC82, VC84, 
VU88, HV91]. Linearized potential analyses are computationally very ef- 
ficient, and for aeroelastic applications, the small disturbance assumption 
is not a severe restriction. However, because of the inherent assumption of 
irrotational and isentropic flow, the linearized potential analyses are only 
valid for subsonic and low transonic speeds and flows without inlet swirl. 
To model the more complex flow fields found in transonic compressors 
and turbines, and as computer speeds and storage have improved, several 
time-accurate time-marching Euler [Hod84, FLM87, Gil88, HSR91] and 
Navier-Stokes [Rai89a, Rai89b, GH91] codes have been developed. These 
codes, unlike the linearized potential codes, account for both steady and 
unsteady vorticity and entropy. However, despite the recent advances in 
computer hardware, the time-marching approach is still computationally 
too expensive for routine design use. 

Recently, several investigators have begun to develop linearized Euler 
analyses [NS76, HC89, HC91, LG91]. As in the linearized potential tech- 
nique, the unsteady flow is assumed to be composed of a nonlinear mean 
flow plus a harmonically varying small perturbation flow. The unsteady 
perturbation flow is governed by the linearized Euler equations. Using this 
model of the flow accounts for entropy and vorticity generation and shocks, 
and permits analysis of unsteady flows about rotational nonisentropic mean 
flows. Furthermore, the solution of the linearized Euler equations, while not 
as efficient as solution of the linearized potential equations, is one to two 
orders of magnitude more efficient than computing the full nonlinear flow 
field while still retaining the essential physics of unsteady flows associated 
with forced response and the onset of flutter. 

In this paper, a recently developed Euler analysis is presented. The 
linearized Euler equations are discretized and solved using a conserva- 
tive pseudo-time time marching Lax-Wendroff scheme [NS76, Ni82, Dan87, 
HC91]. Two important features of the present analysis are the use of a de- 
forming computational grid for flutter calculations, and shock capturing for 
transonic flow calculations. For flutter calculations, the use of a deforming 
computational grid eliminates large error producing gradient terms that 
would otherwise appear in the upwash boundary conditions and in the 
evaluation of the unsteady pressure on the airfoil surfaces. For transonic 
flow calculations, shock capturing is used to model the pressure impulse 
resulting from small harmonic motion of the shock. Shock capturing is 
considerably easier to implement than shock fitting. However, to predict 
correctly the unsteady loads due to the shock impulse, one must use a dis- 
cretization scheme for the unsteady linearized flow solver that is a faithful 
linearization of the nonlinear flow solver [LG91]. Representative computa- 
tional results are presented that demonstrate the capability of the method 
to predict unsteady subsonic and transonic flows. 
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1.2 Theory 


1.2.1 The Linearized Euler Equations 

For compressible flows in which viscous forces may be neglected, the un- 
steady flow is governed by the Euler equations representing the conservation 
of mass, momentum, and energy, i.e, 


dV df dG 

di dx dy 


= 0 


( 1 . 1 ) 


Here, U is the vector of conservation variables and F and G are the so- 
called flux vectors. These are given by 
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Where p is the static density, u and v are the x- and ^components of 
the velocity, p is the pressure, e is the internal energy, and ho is the total 
enthalpy. 

Most previous linearized analyses have used computational grids fixed 
in space. To apply the upwash boundary condition at the airfoil surface, 
the boundary conditions must be extrapolated from the instantaneous po- 
sition of the airfoil to the mean airfoil position, that is, to the boundary 
of the fixed grid. The additional extrapolation terms appearing in the up- 
wash boundary condition contain velocity gradient terms that can be large 
and difficult to evaluate accurately, especially near the leading and trailing 
edges of the airfoils. Whitehead and Grant [WG81] overcame this difficulty 
for rigid-body airfoil motion by using a transformed velocity potential that 
can be viewed as equivalent to using a computational grid that undergoes 
rigid-body motion. Hall [Hal92] has recently developed a linearized poten- 
tial analysis based on a variational principle that includes the effects of a 
grid undergoing arbitrary deformations thus allowing both rigid-body and 
flexible blade motion. 

In the present linearized Euler analysis, we also make use of a deform- 
ing grid. The physical coordinates x,y,t are related to the computational 
coordinates £ , 77 , r by the transformation 



= £ + /(£>»/>r) 

(1.2) 

!/(£)*?> r ) 

= n + g(Z, n, *0 
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<(£> *},t) 
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The grid motion functions / and g are chosen so that the motion of the 
grid conforms to the motion of the airfoils. A typical example of unsteady 
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FIGURE 1.1. Left: Grid in computational coordinate system, r/. Right: Grid 
in physical coordinate system, x,y, for the case of a cascade of airfoils plunging 
with an interblade phase angle, of 90 deg. The cascade shown is Standard 
Configuration No. 10. 


grid motion is shown in Fig. 1.1. Shown is a cascade of airfoils designated 
Standard Configuration No. 10 [Fra91] vibrating in pitch with an interblade 
phase angle, < 7 , of 90 deg. In the computational coordinate system, the 
airfoils and the grid appear to be stationary; in the physical domain, the 
airfoils and grid deform unsteadily. 

Integrating the differential form of the Euler equations [Eq. (1-1)] over a 
control volume D that is fixed in the computational coordinate system - but 
deforming in the physical coordinate system - and applying the divergence 
theorem gives the Euler equations in integral form, i.e., 


(1.5) 


JJ J^Udxdt/H- (f dy - G dx^j — 0 

With some manipulation, Eq. (1.5) can be expressed as 

TrJL^+L [(*-*£) %- I] d " 

Equation (1.6) is valid for both small and large disturbance unsteady 
flows. Large amplitude flows, which exhibit nonlinear behavior, could be 
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computed by time marching Eq. (1.6). However, such computations are 
very expensive. Fortunately, for many flows of interest in aeroelastic ap- 
plications, the unsteadiness in the flow is temporally harmonic and small 
compared with the mean flow. Therefore, to simplify the problem and also 
reduce the required computational time, we assume that the flow is decom- 
posed into a zeroth-order mean or steady flow plus a forst-order unsteady 
perturbation flow. Hence, 

U(f,> 7 ,r) =s U(£,tj) + u(£,??)e' a ' T (1.7) 

Here, U is the vector of conservation variables representing the mean flow 
field, and u is the vector of small perturbation harmonic amplitudes of the 
conservation variables. The motion of the airfoils, and hence the grid, is 
assumed to be small and harmonic so that Eqs. (1.2)-(1.4) are replaced by 



= z+mv)^ 

(1.8) 
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= r 
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where now / and g represent the small perturbation harmonic amplitudes 
of the grid motion. 

Substitution of the perturbation assumptions into the nonlinear Euler 
equations, Eq. (1.6), and collection of terms of zeroth and first-order re- 
sults in the mean flow Euler equations and the linearized Euler equations 
respectively. The mean flow Euler equations in integral form are given by 



— G d£) = 0 


(Ml) 


The linearized Euler equations are given by 

> u II u 1,5 + £> (I? u ' lu" *) = 

— ju f f U ( df dr) + d£ dg) + ju> (/U dr) — g\J d() 

JJd Jsd 
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The Jacobian matrices dF/dV and dG/dU appearing in Eq. (1.12) are 
evaluated using the mean flow solution. For example, for an ideal gas with 
constant specific heats, one finds that the Jacobian dF/dU is given by 

0 1 0 0' 

^U 2 + ^-V 7 -(7-3 )U -(7-l)V 7-1 

—UV v u 0 
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where V ^ = U 2 + V 2 . 

The mean flow [described by (1.11)] is independent of the unsteady per- 
turbation flow. The linearized Euler equations, Eq. (1.12), are linear, in- 
homogeneous, variable coefficient equations. The variable coefficients, i.e., 
the Jacobians, are functions of the mean flow field. The inhomogeneous 
portion of Eq. (M2), that is, the right-hand side, is a function of the mean 
flow and the prescribed grid motion. 

1.2.2 Near-Field Boundary Conditions 

Periodicity 

Because the unsteady flow is governed by linear equations, the response 
of the cascade to blade motion or a gust may be found by decomposing 
the disturbance into a sum of travelling wave modes each with a different 
frequency, w, and interblade phase angle, <r . The total response of the cas- 
cade to the sum of modes is equal to the sum of the responses to each of 
the individual modes. Without loss of generality, we consider here a sin- 
gle travelling wave mode. The complex periodicity condition upstream and 
down stream of the cascade is then 

u«,»7 + G) = u(^v)c j> ( 114 ) 

where G is the blade- to-blade gap. This periodicity condition allows the 
computational domain to be reduced to a single blade passage. 

Flow Tangency 

The upwash boundary condition on the airfoil surface is 

v n = + V • r— ^ * n (Fib) 

where V and v are the mean perturbation flow velocities, n and r are unit 
vectors normal to and tangent to the mean airfoil surface, s is the arc length 
along the airfoil surface, and r is the perturbation displacement of the 
airfoil surface. The two terms on the right-hand side represent the upwash 
due to the local translation and rotation of the airfoil surface. Note that 
because we use a deforming computational grid, no extrapolation terms are 
required. 

1.2.3 Numerical Integration 

To integrate numerically the linearized Euler equations, we use the pseudo- 
time time marching approach suggested by Ni and Sisto [NS76]. The un- 
steady perturbation amplitude u is assumed to vary with time so that 
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Eq. (1.7) is replaced by 

U(£, T,, t) = Ufc, ft) + U«, ft, r)e> UT (1 .16) 

Substitution of Eq. (1.16) into Eq.(1.6) and collection of terms of equal 
order results in the pseudo time-dependent linearized Euler equations, 

irJL^’-Jb^L <U7) 

where b is the same as the right-hand side of the original (time-invariant) 
linearized Euler equation, Eq. (1.12). Note that the first term on the left- 
hand side is the only term that does not appear in the original linearized 
Euler equation. This pseudo- time dependent term is driven to zero by time 
marching the equation to steady state. 

Note the similarity of the linearized Euler equations, Eq. (117), to the 
nonlinear time-dependent Euler equations, Eq. (1.5). One advantage of the 
pseudo-time time-marching approach is that several nonlinear Euler solvers 
have been developed that can be applied with little modification to^the lin- 
earized Euler equations. For example, the conservation variables U in the 
nonlinear Euler solver are replaced by the perturbation conservation vari- 
ables u in the linearized Euler solver. The nonlinear flux vectors F and G 
are replaced by the perturbation flux vectors (9F/3U)u and (dG/d U)u. 
Furthermore, since only the steady state solution of the perturbation solu- 
tion u is desired, acceleration techniques such as multi-grid and local time 
stepping may be used to reduce dramatically the required computational 
time. 

For subsonic flow calculations, any consistent discretization of the lin- 
earized Euler equations is sufficient to produce consistent unsteady flow 
solutions. For transonic flows with shock capturing, however, the discretized 
perturbation equations should be a faithful linearization of the discretized 
unsteady nonlinear equations used to compute the steady flow [LG91]. In 
the present study, we use Ni’s conservative finite volume implementation of 
the Lax-Wendroff method [Ni82, Dan87] to compute both the steady and 
unsteady flow fields. Local time stepping and Ni’s multiple-grid accelera- 
tion technique are used to speed convergence in both the steady nonlinear 
Euler and unsteady linearized Euler codes. For subsonic flows, a fourth- 
difTerence smoothing operator is used to eliminate sawtooth modes. For 
transonic flows, a combination of fourth- difference and second difference 
smoothing is used to capture shocks. 


1.2.4 Far- Field Boundary Conditions 

Because the computational domain extends a finite distance in the axial 
direction (see Fig. 1.1), far-field boundary conditions are required at the 
far-field boundaries to prevent spurious reflections of outgoing waves. Pre- 
vious investigators have found analytical expressions for the behavior of 
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the linearized equations in the far field, and have matched these far-field 
solutions to the near-field numerical solutions at the upstream and down- 
stream far-field boundaries. These analytical solutions were found for the 
special case of uniform mean flow in the far field. For transonic Euler flows, 
however, the mean flow in the far field is in general nonuniform. For ex- 
ample, the presence of a shock or multiple shocks in the blade passage will 
produce nonuniform total pressure in the downstream region. No general 
analytical expression has been developed for this situation. 

To overcome this difficulty, the authors have developed an alternative 
treatment the far-field boundary conditions. Instead of finding the analyt- 
ical eigenmodes of the linearized Euler equations in the far-field, we find 
the eigenmodes of the discretized linearized Euler equations. The numer- 
ical eigenmodes are then matched to the numerical near-field solution at 
the far-field boundaries. This procedure, although computationally less ef- 
ficient, produces numerically exact nonreflecting boundary conditions. Ref- 
erence [HC91] describes the implementation of these boundary conditions 
in more detail. 


1.3 Results 

1.3.1 Transonic Channel Flow 

The first test case to be presented is of transonic flow through a diverging 
channel. This case is presented to demonstrate the ability of the linearized 
Euler method to model accurately shock motion using shock capturing. 
The channel considered here has a height, A ) given by 

A(x) = j4j n i et jl. 10313+ .10313 tanh [ 5 (* “ ^)]} > 0 < * < 1 ( 118 ) 

So that the results can be compared to a one-dimensional shock-fitting 
theory, A; n ] et is taken to be small compared with the channel length 
(Ainiet = 0 01)- The inflow total pressure, Pt, total density, pr, and flow 
velocity, U } are 1.0, 1.364, and 1.0, respectively. The back pressure, P ? xit is 
0.7422. Shown in Fig. 1.2 is the Mach number and pressure distribution as 
computed by the present nonlinear steady Euler solver using a 129 x 5 node 
computational grid. Also shown for comparison is the solution found using 
a steady quasi-one-dimensional, shock-fitting, Euler solver using 1001 grid 
nodes in the x direction. Note the excellent agreement between the two 
approaches. The only noticeable differences occur at the shock where the 
present nonlinear Euler solver smears the shock over about five grid nodes. 

Next, we consider the case of a quasi-steady perturbation in the back 
pressure. The perturbation solution was computed using three different 
methods. We first computed the solution using the present linearized solver; 
then using a quasi one-dimensional, shock-fitting, linearized Euler solver; 
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FIGURE 1.2. Steady transonic flow in a diverging channel. 


and finally, by differencing the solutions found using the present nonlinear 
steady Euler solver at two slightly different back pressures. The results 
are shown in Fig. 1.3. The three solutions are seen to be in very good 
agreement away from the shock. At the shock, the two shock-capturing 
techniques produce an impulse of pressure. This impulse represents the 
load exerted on the channel walls due to the displacement of the shock. 
Note that the present linearized Euler solution is nearly identical to the 
differenced nonlinear solutions. 

We now consider the case of an unsteady variation in the back pres- 
sure. Shown in Fig. 1.4 are the computed real and imaginary parts of the 
unsteady pressure distribution due to a unit perturbation in the back pres- 
sure with an excitation frequency, w, of 1.0. The results of the present lin- 
earized Euler code agree well with the quasi one-dimensional, shock-fitting, 
linearized Euler solver. 

To determine whether the present linearized Euler solver correctly pre- 
dicts the unsteady loads induced by the motion of the shock, we integrated 
the pressure along the lower wall to obtain the “wall force.” These results 
are tabulated in Table 1.1 for several different excitation frequencies. Also 
given in Table 1.1 is the wall force computed using the shock-fitting code. 
For our purposes, the grid is sufficiently fine that the shock-fit solution 
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FIGURE 1.3. Perturbation pressure in a diverging channel due to a steady per- 
turbation in back pressure. Bottom figure shows detail of shock impulse. 


may be considered exact. For low frequencies, the two codes are in almost 
perfect agreement suggesting that the shock impulse found using shock cap- 
turing is properly modelled. For the higher frequency cases, the agreement 
is still good, but there is a slight error in the phase of the wall force. We 
believe that these differences are due to dispersion errors in regions away 
from the shock rather than a limitation in shock capturing at high frequen- 
cies. Finally, For the uj = 1 case, we deliberately introduced two types of 
inconsistency into the linearized Euler solver to determine their influences 
on the solutions. In the first case, we used a different level of smoothing in 
the linearized unsteady Euler code than was used in the nonlinear steady 
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FIGURE 1.4. Unsteady pressure in a diverging channel due to an unsteady per- 
turbation in back pressure. The reduced frequency, u>, based on channel length 
and inflow velocity is 1.0. 

Euler code. In the second case, global time stepping was used in the steady 
Euler code and local time stepping was used in the linearized unsteady 
code. In both cases, the inconsistency between steady and unsteady codes 
is seen to increase the error in the unsteady wall force. The errors are much 
more pronounced, however, for the latter case. These results are consistent 
with the claim by Lindquist and Giles that the linearized code should be a 
faithful linearization of the nonlinear code [LG91]. 

1.3.2 Subsonic Flow in a Compressor 

Having demonstrated the accuracy of the method for transonic channel 
flow, we next consider unsteady flows in a compressor geometry to demon- 
strate the ability of the method to accurately compute the unsteady flow 
about a cascade of loaded airfoils. The cascade considered is the newly des- 
ignated Standard Configuration No. 10 [Fra91]. The airfoils have a circular 
arc camber distribution with a maximum height of 5 percent of the chord 
with a slightly modified NACA 0006 thickness distribution. The stagger 
angle, 6, is 45 deg and the gap-to-chord ratio, G, is 1 . 0 . For the subsonic 
case considered here, the steady inflow angle, f!_oo> i 8 55 deg and the in- 
flow Mach number, M_ , is 0.7. Under these conditions, the flow is entirely 
subsonic with a maximum Mach number on the suction surface of about 
0.92. 

Figure 1.5 shows the computed unsteady pressure distribution for the 
case where the airfoils pitch about points near their midchords with an 
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TABLE 1.1. Predicted pressure loads in a transonic diverging channel due to an 
unsteady perturbation in back pressure. 


Frequency, u> Scheme Wall Force 

0lT~ ID Shock Fitting 1.0305 + 0.0000;” 

Nonlinear Euler° 1.0326 + 0.0000; 

Present Method 1.0320 + 0.000 0; 

T5 ID Shock Fitting 0.6866 + 0.6115; 

Present Method 0.6863 + 0.6100; 

LO ID Shock Fitting 0.1249 — 0.6267;” 

Present Method 0.1293 — 0.6222; 

Present Method** 0.1315 — 0.6184; 

Present Method c 0.1778 — 0.6133; 

l0 ID Shock Fitting -0.0807 - 0.180lJ 

Present Method —0.0750 — 0.1831; 


Shock Force 
0.5133 + 0.0000; 


0.1796 — 0.4426; 
-0.2543 - 0.2797J 


-0.1182 + 0.1667; 


° Results from the steady analysis were found for two slightly different back pressures. The 
two solutions were then differenced and normalized by AP ex jj; . . - 

^Smoothing coefficient in steady analysis is 0.15; smoothing coefficient in unsteady analysis 

18 0 ^Giobal time stepping used in steady analysis; local time stepping used in unsteady anal- 


ysis. 


interblade phase angle, <r, of 90 deg and a reduced frequency, w, of 0.5 
(Standard Configuration No. 10, Case 2). Here, the pressure is nondimen- 
sionalized by the upstream steady flow quantity pV % - To assess the accuracy 
of the analysis, solutions were computed on both a 65 x 17 and a 129 x 33 
node grid. The good agreement between the coarse and fine grid solutions 
indicates that the fine grid solution is nearly grid converged. Also shown for 
comparison is the solution found using a deforming grid linearized potential 
method [Hal92]. Note that the solutions are well behaved in the neighbor- 
hood of the leading and trailing edges. The use of a deforming grid greatly 
improves the accuracy of the solution in these regions. 

The linearized Euler solver is very efficient. The unsteady solution com- 
puted on the 65 x 17 node grid required 565 seconds of computer time on an 
IBM RISC System/6000 model 320 computer. For comparison, the steady 
solution required 110 seconds. 

Having computed the unsteady pressure distribution on the airfoil sur- 
face, one can integrate to obtain the pitching moment. Shown in Figure 1.6 
is the imaginary part of the pitching moment computed for a range of 
interblade phase angles, -180 < a < 180, at a reduced frequency, w, of 
1.0. Also shown is the pitching moment computed using Hall’s deforming 
grid linearized potential method. The peaks in the solution are acoustic 
resonances. We were unable to obtain converged solutions near one of the 
upstream acoustic resonance points (<r = 117.1 deg). Note that there is gen- 
erally good agreement between the linearized Euler and linearized potential 
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Chordal Distance, x’/c 

FIGURE 1.5. Unsteady pressure on surface of subsonic compressor blade due to 
pitching motion (Standard Configuration No. 10, Case 2). The reduced frequency, 
u7, is 0.5, and the interblade phase angle, <r, is 90 deg. 


solution. One noticeable difference is seen at c = —20 deg. This point is 
very near one of the upstream acoustic resonance points (< t = —26.93). 
Nevertheless, the good agreement elsewhere in the superresonant regions 
demonstrates the effectiveness of the numerical far-field boundary condi- 
tions. 





14 1. Calculation of Unsteady Linearized Euler Flows in Cascades 



' 180 . *- 120 . - 60 . 0 . 60 . 120 . 180 , 
Interblade Phase Angle (deg) 


FIGURE 1.6. Imaginary pitching moment acting on subsonic compressor blades 
due to pitching motion (Standard Configuration No. 10). The reduced frequency, 
u7, is 1.0. 


1.3.3 Transonic Flow in a Compressor 

The final case to be considered here is flow through a transonic compressor. 
The cascade is again Standard Configuration No. 10. However, the flow con- 
ditions are such that there is a supersonic patch on the suction surface of the 
airfoil. The mean inflow angle, , is 58 deg and the inflow Mach number, 
M-oo , is 0.8. Shown in Fig. 1,7 is the computed isentropic Mach number 
distribution on the surface of the airfoil. Also shown is the Mach number 
distribution computed using a full potential technique [VC84, VU88]. Note 
that for the solution computed with the present nonlinear Euler solver, the 
shock is smeared over about five grid points. 

Having computed the steady flow solution, we consider the case of the 
airfoils vibrating in pitch with an interblade phase angle, <x, of 90 deg 
with a reduced frequency, uJ, of 1.0 (Standard Configuration No. 10, Case 
20). Figure 1.8 shows the computed unsteady pressure distribution on the 
airfoil surface. The shock impulse representing the unsteady load due to 
unsteady shock motion is clearly visible on the suction surface. Also shown 
is the pressure distribution computed using Verdon’s linearized potential 
code, LINSUB code [VU88], The linearized Euler and linearized potential 
solutions agree quite well on the pressure surface. However, the agreement 
is not as good on the suction surface, particularly ahead of the shock. This 
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Chordal Distance, x’/c 


FIGURE 1.7. Steady pressure distribution on transonic compressor blade (Stan- 
dard Configuration No. 10). The inflow Mach number, M , is 0.8, the gap-to-chord 
ratio, G, is 1.0, the stagger angle, 0, is 45 deg., and the inflow angle, fl-oo, is 58 
deg. 


difference could arise from neglecting entropy and vorticity generation in 
the steady and unsteady potential analyses. 


1.4 Summary 

A two-dimensional linearized Euler solver has been developed that is capa- 
ble of computing unsteady small disturbance flows in cascades. Two impor- 
tant features of the analysis are the use of a deforming grid to improve the 
accuracy of flutter calculations, and the use of shock capturing to model 
shock motion. Using the pseudo-time time-marching technique, one can 
solve the linearized equations with any of a number of well known schemes 
originally developed to solve the nonlinear Euler equations. Furthermore, 
since only the steady state harmonic perturbation solution is desired, con- 
vergence acceleration techniques may be used to speed convergence. Re- 
sults presented in the paper clearly demonstrate the ability of the analysis 
to predict unsteady subsonic and transonic flows accurately and efficiently. 
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FIGURE 1.8. Unsteady pressure on surface of transonic compressor blade due 
to pitching motion (Standard Configuration No. 10, Case 20). The reduced fre- 
quency, w, is 1.0, and the interblade phase angle, <r, is 90 deg. 
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